Preprint typeset in JHEP style - HYPER VERSION 



UFIFT-HEP-05-05 
MCTP-05-42 



Contrasting Supersymmetry and Universal Extra 
Dimensions at the CLIC Multi-TeV e + e~~ Collider 



\ Marco Battaglia 

o . 

Dept. of Physics, University of California, Berkeley, CA 94720, USA 
<N ' and 



(N 

> 



Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA 



E-mail: MBattaglia@lbl.gov 



Asesh Krishna Datta 

MCTP, University of Michigan, Ann Arbor, MI 48109-1120, USA 



E-mail: 



aseshOumich . edu 



Tj- . Albert De Roeck 



CERN, Geneva, Switzerland 



E-mail: Albert . de . Roeck@cern . ch 



o 

O 
in 

O ' Kyoungchul Kong 

Physics Dept., University of Florida, Gainesville, FL 32611, USA 
E-mail: kong@phys.ufl.edu 

Oh! 

D . Konstantin T. Matchev 



Physics Dept., University of Florida, Gainesville, FL 32611, USA 
E-mail: natchev@phys . uf 1 . edu| 



X 

Abstract: Universal extra dimensions and supersymmetry have rather similar experi- 
mental signatures at hadron colliders. The proper interpretation of an LHC discovery in 
either case may therefore require further data from a lepton collider. In this paper we 
identify methods for discriminating between the two scenarios at the linear collider. We 
study the processes of Kaluza-Klein muon pair production in universal extra dimensions in 
parallel to smuon pair production in supersymmetry, accounting for the effects of detector 
resolution, beam-beam interactions and accelerator induced backgrounds. We find that the 
angular distributions of the final state muons, the energy spectrum of the radiative return 
photon and the total cross-section measurement are powerful discriminators between the 
two models. Accurate determination of the particle masses can be obtained both by a 
study of the momentum spectrum of the final state leptons and by a scan of the particle 
pair production thresholds. We also calculate the production rates of various Kaluza-Klein 
particles and discuss the associated signatures. 
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1. Introduction 

Supersymmetry (SUSY) and Extra Dimensions offer two different paths to a theory of new 
physics beyond the Standard Model (SM). They both address the hierarchy problem, play 
a role in a more fundamental theory aimed at unifying the SM with gravity, and offer 
a candidate particle for dark matter, compatible with present cosmology data. If either 
supersymmetry or extra dimensions exist at the TeV scale, signals of new physics should be 
found by the ATLAS and CMS experiments at the Large Hadron Collider (LHC) at CERN. 
However, as we discuss below, the proper interpretation of such discoveries, namely the 
correct identification of the nature of the new physics signals, may not be straightforward 
at the LHC and may require the complementary data from an e + e~ collider. In particular, 
a multi-TeV collider, such as CLIC [1,2], would ensure a sensitivity over a broad mass 
range. 

A particularly interesting scenario of TeV-size extra dimensions is offered by the so 
called Universal Extra Dimensions (UED), originally proposed in [3], where all SM particles 
are allowed to freely propagate into the bulk. The case of UED bears interesting analogies 
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to supersymmetry 1 , and sometimes has been referred to as "bosonic supersymmetry" [4]. 
In principle, disentangling UED and supersymmetry appears highly non-trivial at hadron 
colliders [4-6]. For each SM particle, both models predict the existence of a partner (or 
partners) with identical interactions. Unfortunately, the masses of these new particles 
are model-dependent and cannot be used to unambiguously discriminate between the two 
theories, although a degenerate spectrum might suggest UED while a split spectrum would 
hint 2 towards SUSY. One would therefore like to have an experimental verification which 
relies on the fundamental distinctions between the two models. 

One of the characteristic features of UED is the presence of a whole tower of Kaluza- 
Klein (KK) partners, labelled by their KK level n. In contrast, N = 1 supersymmetry 
predicts a single superpartner for each SM particle. One might therefore hope to be able 
to discover the higher KK modes of UED and, having observed a repetition of the SM 
particle content at each KK level, prove the existence of extra dimensions. However, there 
are two significant challenges along this route. First, the masses of the higher KK modes 
are (roughly) integer multiples of the masses of the n = 1 KK partners, and as a result their 
production cross-sections are kinematically suppressed. Second, they predominantly decay 
to n = 1 KK modes, thus simply contributing a small amount to the inclusive production 
of n = 1 KK particles. Just as in the case of SUSY, because of the unknown momentum 
carried away by the dark matter candidate at the end of the decay chain, one is unable 
to reconstruct individual KK resonances. The only exceptions are the level 2 KK gauge 
bosons, which may appear as high mass dijet or dilepton resonances, when they decay 
directly to SM fermions through loop suppressed couplings [4,10-12]. Unfortunately the 
reach of the LHC is rather limited in this case [13], due to the smallness of these loop- 
suppressed branching fractions. In addition, even if they are discovered, such resonances 
can be misinterpreted as ordinary Z' gauge bosons in extended supersymmetric models [13]. 

The second fundamental distinction between SUSY and the minimal version o/UED is 
that supersymmetric extensions of the SM necessarily have extended Higgs sectors leading 
to additional Higgs and higgsino states in the spectrum. The higgsinos of supersymmetry 
are in one-to-one correspondence with the n = 1 KK partners of the Goldstone bosons and 
cannot be used for discrimination. The additional Higgs bosons of SUSY are not a robust 
discriminator either. On the one hand, they may simply escape detection at the LHC: the 
SUSY parameter space has large portions where the LHC discovers a single (SM-like) Higgs 
boson and misses the rest of the Higgs spectrum. Alternatively, one may consider UED 
models with an extended Higgs sector which would mimic the SUSY Higgs phenomenology. 

The third and most fundamental distinction between UED and supersymmetry is re- 
flected in the properties of the individual particles: the KK partners have identical spin 

1 More precisely, the phenomenology of the first level (n = 1) of Kaluza-Klein modes in UED is very 
similar to that of N = 1 supersymmetric models with somewhat degenerate superpartner spectrum and 
stable lightest supersymmetric particle (LSP). In what follows we shall use the term "supersymmetry" in 
this somewhat narrower context. 

2 Notice that the recently proposed little Higgs models with T-parity [7-9] are reminiscent of UED, and 
their spectrum does not have to be degenerate, so they may still be confused with supersymmetry. Under 
those circumstances, the methods for discrimination discussed in this paper would still apply. 
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quantum numbers as their SM counterparts, while the spins of the superpartners differ 
by 1/2 unit. However, spin determinations also appear to be difficult at the LHC (or at 
hadron colliders in general), where the center of mass energy in each event is unknown. In 
addition, the momenta of the two dark matter candidates in the event are also unknown. 
This prevents the reconstruction of any rest frame angular decay distributions, or the di- 
rections of the two particles at the top of the decay chains. Recently it has been suggested 
that a charge asymmetry in the lepton-jet invariant mass distributions from a particular 
cascade, can be used to discriminate SUSY from the case of pure phase space decays [14]. 
However, the possibility of discriminating SUSY and UED by this method still needs to be 
demonstrated. 

At the LHC, strong processes dominate the production of both KK modes and super- 
partners. The resulting signatures involve relatively soft jets, leptons and missing transverse 
energy [4]. The study in [4] shows that, with 100 fb _1 of integrated luminosity, the LHC 
experiments will be able to cover all of the cosmologically preferred parameter space in 
UED [15]. However, due to the reasons explained above, the discrimination between UED 
and supersymmetry may be difficult at the LHC. 

In this paper, we focus on the minimal UED model, described in some detail in Sec- 
tion ||, and on the Minimal Supersymmetric extension of the Standard Model (MSSM) 
to study how the two models can be discriminated at a linear collider tunable over a 
center-of-mass energy range of 1 TeV to 3 TeV. We assume that the LHC will have already 
observed signals of new physics consistent with either n = 1 KK modes in UED or sparticle 
production in supersymmetry. We concentrate on the relevant questions to be addressed 
by a post-LHC facility. In section |3| we describe the event simulation and reconstruction 
adopted in this analysis, while in section || we contrast UED with supersymmetry using 
the example of y&pTTfj final state. In section || we extend our discussion to other possible 
final states. Section || has our conclusions. 

2. The minimal UED model 

In its simplest incarnation, the UED model has all the SM particles propagating in a single 
extra dimension of size R, which is compactified on an S1/Z2 orbifold. More complicated 
versions have also been proposed, motivated by ideas about electroweak symmetry break- 
ing [16], neutrino masses [17,18], proton stability [19] or the number of generations [20]. 
A peculiar feature of UED is the conservation of Kaluza-Klein number at tree level, which 
is a simple consequence of momentum conservation along the extra dimension. However, 
bulk and brane radiative effects [10-12] break KK number down to a discrete conserved 
quantity, the so called KK parity, (—1)™, where n is the KK level. KK parity ensures that 
the lightest KK partners - those at level one - are always pair-produced in collider exper- 
iments, similar to the case of supersymmetry models with conserved i?-parity. KK parity 
conservation also implies that the contributions to various precisely measured low-energy 
observables [21-28] only arise at loop level and are small. As a result, the limits on the scale 
of the extra dimension, from precision electro- weak data, are rather weak, constraining R^ 1 
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to be larger than approximately 250 GeV [25]. An attractive feature of UED is the pres- 
ence of a stable massive particle which can be a cold dark matter candidate [12,15,29,30]. 
The lightest KK partner (LKP) at level one has negative KK parity and can be stable on 
cosmological scales. In the minimal model, the LKP is the KK partner of the hypercharge 
gauge boson [12] and its relic density is generically in the desired range [15]. Kaluza-Klein 
dark matter offers excellent prospects for direct or indirect detection [30-36]. Once the 
radiative corrections to the Kaluza-Klein masses are properly taken into account [12], the 
collider phenomenology of the minimal UED model 3 exhibits striking similarities to super- 
symmetry [4, 38] and represents an interesting and well motivated counterexample which 
can "fake" supersymmetry signals at the LHC. 

In the minimal UED model, the bulk interactions of the KK modes readily follow from 
the SM Lagrangian and contain no unknown parameters other than the mass, Mh, of the 
SM Higgs boson. In contrast, the boundary interactions, which are localized on the orbifold 
fixed points, are in principle arbitrary, and thus correspond to new free parameters in the 
theory. They are in fact renormalised by bulk interactions, and are scale dependent [10]. 
Therefore, we need an ansatz for their values at a particular scale. Of course, the UED 
model should be treated only as an effective theory which is valid up to some high scale 
A, at which it is matched to some more fundamental theory. In the minimal UED model 
the boundary terms are assumed to vanish at the cutoff scale A, and are subsequently 
generated through RGE evolution to lower scales. Thus the minimal UED model has only 
two input parameters: the size of the extra dimension, R, and the cutoff scale, A. The 
number of KK levels present in the effective theory is therefore AR and may vary between 
a few and ~ 40, where the upper limit corresponds to values of A leading to a breakdown 
of perturbativity below the A scale. 

3. Event simulation and data analysis 

In order to study the discrimination of UED signals from supersymmetry, we have im- 
plemented the relevant features of the minimal UED model in the CompHEP event gen- 
erator [39]. The MSSM is already available in CompHEP since version 41.10. All n = 1 
KK modes are incorporated as new particles, with the proper interactions and one-loop 
corrected masses [12]. The widths can then be readily calculated with CompHEP on a case 
by case basis and added to the particle table. Similar to the SM case, the neutral gauge 
bosons at level 1, Z\ and 71, are mixtures of the KK modes of the hypercharge gauge boson 
and the neutral SU{2)w gauge boson. However, it was shown in [4] that the radiatively 
corrected Weinberg angle at level 1 and higher is very small. For example, 71, which is the 
LKP in the minimal UED model, is mostly the KK mode of the hypercharge gauge boson. 
For simplicity, in the code we neglect neutral gauge boson mixing for n > 1. 

In the next section we concentrate on the pair production of level 1 KK muons e + e~ — > 
and compare it to the analogous process of smuon pair production in supersymmetry: 
e+e~ — > jj- + fi~. In UED there are two n = 1 KK muon Dirac fermions: an SU(2)\y doublet 

3 For alternative possibilities, see [37]. 
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(a) 



(b) 



Figure 1: The dominant Feynman diagrams for KK muon production e + e — ► /if \i x in Universal 
Extra Dimensions. The black dot represents a KK-number violating boundary interaction [12]. 



Hi and an SU(2)\y singlet fx x , both of which contribute in Eq. (|3.l| ) below (see also Fig. 
In complete analogy, in supersymmetry, there are two smuon eigenstates, /2l and £lr, both 
of which contribute in Eq. (|3.2| ). The dominant diagrams in that case are shown in Fig. |2|. 




Figure 2: 



In principle, there are also diagrams mediated by 
7„,Z n for n = 4,6,... but they are doubly suppressed 
- by the KK-number violating interaction at both ver- 
tices and the KK mass in the propagator - and here 
can be safely neglected. However, 72 and Z2 exchange 
(Fig. |]b) may lead to resonant production and signifi- 
cant enhancement of the cross-section, as well as inter- 



esting phenomenology as discussed below in Section 4.5 



The dominant Feyn- 
man diagrams for smuon production 
e + e~ — > in supersymmetry. 



We have implemented the level 2 neutral gauge bosons 
72,^2 with their widths, including both KK-number 
preserving and the KK-number violating decays as in 
Ref. [4]. We consider the final state consisting of two opposite sign muons and missing 
energy. It may arise either from KK muon production in UED 



e e -> n{p x -> h [i 71 71 , 
with 71 being the LKP, or from smuon pair production in supersymmetry: 



(3.1) 



e e 



(3.2) 



where x± is the lightest supersymmetric particle. We reconstruct the muon energy spectrum 
and the muon production polar angle, aiming at small background from SM processes with 
minimal biases due to detector effects and selection criteria. The goal is to disentangle KK 
particle production ( |3.l| ) in UED from smuon pair production (3^) in supersymmetry. We 
also determine the masses of the produced particles and test the model predictions for the 
production cross-sections in each case. 

We first fix the UED parameters to R~ l = 500 GeV, AR = 20, leading to the spectrum 
given in Table [l]. 
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Particle 


Mass 


Mi 


515.0 GeV 


Mi 


505.4 GeV 


7i 


500.9 GeV 



Table 1: Masses of the KK excitations for the parameters R 1 = 500 GeV and AR = 20 used in 
the analysis. 

The ISR-corrected signal cross section in UED for the selected final state p + LT 7171 is 
14.4 fb at sj~s = 3 TeV. Events have been generated with CompHEP and then reconstructed 
using a fast simulation based on parametrized response for a realistic detector at CLIC. 
In particular, the lepton identification efficiency, momentum resolution and polar angle 
coverage are of special relevance to this analysis. We assume that particle tracks will be 
reconstructed through a discrete central tracking system, consisting of concentric layers of 
Si detectors placed in a 4 T solenoidal field. This ensures a momentum resolution 5p/p 2 
= 4.5 xlO -5 GeV -1 . A forward tracking system should provide track reconstruction down 
to ~ 10°. We also account for initial state radiation (ISR) and for beamstrahlung effects 
on the center-of-mass energy. We assume that muons are identified by their penetration in 
the instrumented iron return yoke of the central coil. A 4 T magnetic field sets an energy 
cutoff of ~ 5 GeV for muon tagging. 

The events from the CompHEP generation have been treated with the Pythia 6 . 210 par- 
ton shower [40] and reconstructed with a modified version of the SimDet 4 . program [41]. 
Beamstrahlung has been added to the CompHEP generation. The luminosity spectrum, ob- 
tained by the GuineaPig beam simulation for the standard CLIC beam parameters at 
3 TeV, has been parametrised using a modified Yokoya-Chen approximation [42]: 

This analysis has backgrounds coming from SM /i + li~vv final states, which are mostly 
due to gauge boson pair production W + W~ — > pT v Z®Z® — ► p + p~vv and from 
e + e~ — ► W + W~u e u e , e + e~ — > Z°Z°i/ e i/ e , followed by muonic decays. The background 
total cross section is ~20 fb at y/s = 3 TeV. In addition to its competitive cross section, 
this background has leptons produced preferentially at small polar angles, therefore biasing 
the angular distribution. In order to reduce this background, a suitable event selection has 
been applied. Events have been required to have two muons, missing energy in excess to 
2.5 TeV, transverse energy below 150 GeV and event sphericity larger than 0.05. In order 
to reject the Z°Z° background, events with di-lepton invariant mass compatible with M z o 
have also been discarded. The underlying 77 collisions also produce a potential background 
to this analysis in the form of 77 — > This background has been simulated using the 

CLIC beam simulation and Pythia. Despite its large cross section, it can be completely 
suppressed by a cut on the missing transverse energy E™ tssm 9 > 50 Q e y. Finally, in order 
to remove events with large beamstrahlung, the event sphericity had to be smaller than 
0.35 and the acolinearity smaller than 0.8. These criteria provide a factor ~ 30 background 
suppression, in the kinematical region of interest, while not significantly biasing the lepton 
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momentum distribution. 



4. Comparison of UED and supersymmetry in fi + /j, Ifi 

In order to perform the comparison of UED and MSSM, we adjusted the MSSM parameters 
to get the two smuon masses M^ L and Mf LR and the lightest neutralino mass M^o matching 
exactly those of the two Kaluza-Klein muons M^p and M^s and of the KK photon M 71 for 
the chosen UED parameters. It must be stressed that such small mass splitting between 
the two muon partners is typically rather accidental in supersymmetric scenarios. The 
supersymmetric parameters used are given in Table || 



MSSM Parameter 


Value 


A* 


1000 GeV 


Afi 


502.65 GeV 


M 2 


1005.0 GeV 


M» L 


512.83 GeV 


M^ R 


503.63 GeV 


tan/3 


10 



Table 2: MSSM parameters for the SUSY study point used in the analysis. This choice of soft 
SUSY parameters in CompHEP Jeads to an exact match between the corresponding UED and SUSY 
mass spectra. 

We then simulate both reactions ( [3.1D and ( |3.2| ) with CompHEP and pass the resulting 
events through the same simulation and reconstruction. The ISR-corrected signal cross- 
section in SUSY for the selected final state /J- + /J-~XiXi 1S 2.76 fb at \fs= 3 TeV, which is 
about 5 times smaller than in the UED case. 



4.1 Angular distributions and spin measurements 

In the case of UED, the KK muons are fermions and their angular distribution is given by 

dcos9J UED E^+M^ 

Assuming that at CLIC the KK production takes place well above threshold, the formula 
simplifies to: 



da 



~l + cos 2 6>. (4.2) 

UED 



K d cos i 

As the supersymmetric muon partners are scalars, the corresponding angular distribution 



is 

(4.3) 

dcos8 J SUSY 
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Figure 3: Differential cross-section da jd cos 8^ for UED (blue, top) and supersymmetry (red, 
bottom) as a function of the muon scattering angle 9^. The figure on the left shows the ISR- 
corrected theoretical prediction. The two figures on the right in addition include the effects of event 
selection, beamstrahlung and detector resolution and acceptance. The left (right) panel is for the 
case of UED (supersymmetry). The data points are the combined signal and background events, 
while the yellow-shaded histogram is the signal only. 

Distributions ( |4.2| ) and ( |4.3| ) are sufficiently distinct to discriminate the two cases. 
However, the polar angles of the original KK-muons and smuons are not directly observ- 
able and the production polar angles 9^ of the final state muons are measured instead. But 
as long as the mass differences M Ml — M 7l and — M^o respectively remain small, the 
muon directions are well correlated with those of their parents (see Figure |j|a) . In Fig. [||b 
we show the same comparison after detector simulation and including the SM background. 
The angular distributions are well distinguishable also when accounting for these effects. 
By performing a \ 2 fit to the normalised polar angle distribution, the UED scenario con- 
sidered here could be distinguished from the MSSM, on the sole basis of the distribution 
shape, with 350 fb _1 of data at = 3 TeV. 

4.2 Threshold scans 

At the e + e~ linear collider, the muon excitation masses can be accurately determined 
through an energy scan of the onset of the pair production threshold. This study not only 
determines the masses, but also confirms the particle nature. In fact the cross sections for 
the UED processes rise at threshold oc (3 while in supersymmetry their threshold onset is 
oc /3 3 , where (5 is the particle velocity. 

Since the collision energy can be tuned at properly chosen values, the power rise of the 
cross section can be tested and the masses of the particles involved measured. We have 
studied such threshold scan for the e + e~ — » fJ^^i — > /i + ^~7i7i process at yfs = 1 TeV, 
for the same parameters as in Table |]. We account for the anticipated CLIC centre-of-mass 
energy spread induced both by the energy spread in the CLIC linac and by beam-beam 
effects during collisions. This been obtained from the detailed GuineaPig beam simulation 
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Figure 4: The total e + e~ — > /Ui Mi — ► /i + /j~7i7i cross section a in pb as a function of the center- 
of-mass energy y/s near threshold. Left: the threshold onset with (line, blue) and without (dots) 
beamstrahlung effects. Right: a threshold scan at selected points. The green curve refers to the 
reference UED parameters while for the red (blue) curve the mass of /if (Hi) has been lowered by 
2.5 GeV. The points indicate the expected statistical accuracy for the cross section determination 
at the points of maximum mass sensitivity. Effects of the CLIC luminosity spectrum are included. 

and parametrised using the modified Yokoya-Chen model [42,43]. An optimal scan of 
a particle pair production threshold consists of just two energy points, sharing the total 
integrated luminosity in equal fractions and chosen at energies maximising the sensitivity 
to the particle widths and masses [44]. For the UED model scan we have taken three 
points, one for normalisation and two at the maxima of the mass sensitivity (see Figure Q). 
Inclusion of beamstrahlung effects induces a shift of the positions of these maxima towards 
higher nominal -y/s values [45] . From the estimated sensitivity da jdM and the cross section 
accuracy, the masses of the two UED muon excitations can be determined to ±0.11 GeV 
and ±0.23 GeV for the singlet and the doublet states respectively, with a total luminosity 
of 1 ab _1 shared in three points, when the particle widths can be disregarded. 

4.3 Production cross section determination 

The same analysis can be used to determine the cross section for the process e + e _ — > 
The SM contribution can be determined independently, using anti-tag cuts, and 
subtracted. Since the cross section for the UED process at 3 TeV is about five times 
larger compared to smuon production in supersymmetry, this measurement would reinforce 
the model identification obtained by the spin determination. This can be quantified by 
performing the same x 2 fit to the muon polar production angle discussed above, but now 
including also the total number of selected events. Since the cross section depends on the 
mass of the pair produced particles, we include a systematic uncertainty on the prediction 
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Figure 5: The muon energy spectrum resulting from KK muon production (3.1) in UED (blue, top 
curve) and smuon production ( |3.2|) in supersymmetry (red, bottom curve). The UED and SUSY 
parameters are chosen as in Fig. ^. The plot on the left shows the ISR-corrected distribution, while 
that on the right includes in addition the effects of event selection, beamstrahlung and detector 
resolution and acceptance. The data points are the combined signal and background events, while 
the yellow-shaded histogram is the signal only. 



corresponding to a ±0.05 % mass uncertainty, which is consistent with the results discussed 
below. At CLIC the absolute luminosity should be measurable to O(0.1 %) and the average 
effective collision energy to 0(0.01 %). 

4.4 Muon energy spectrum and mass measurements 

The characteristic end-points of the muon energy spectrum are completely determined by 
the kinematics of the two-body decay and hence they don't depend on the underlying 
framework (SUSY or UED) as long as the masses involved are tuned to be identical. We 
show the ISR-corrected expected distributions for the muon energy spectra at the generator 
level in Fig. ||a, using the same parameters as in Fig. ||. As expected, the shape of the 
distribution in the case of UED coincides with that for MSSM. 

The lower, E m i n , and upper, E max , endpoints of the muon energy spectrum are related 
to the masses of the particles involved in the decay according to the relation: 

1 / M2 o \ 

E max /rmn = ^M fi ( 1 - 1 7 (1 ± 0) (4.4) 

where Mjx and M^o are the smuon and LSP masses and 7 = 1/(1 — /3 2 ) 1 / 2 with (3 = 

Jl - M?/E% eam is the p, boost. In the case of the UED the formula is completely analogous 
with M w replacing and M 71 replacing M^o . 

Due to the splitting between the fii and fiR masses in MSSM and that between the 
/Ltf and nf masses in UED, in Fig. [|a we see the superposition of two box distributions. 
The left, narrower distribution is due to /if pair production in UED pair production in 
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supersymmetry). The underlying, much wider box distribution is due to /if pair production 
in UED (Jxl pair production in supersymmetry). The upper edges are well defined, with 
smearing due to beamstrahlung and, but less importantly, to momentum resolution. The 
lower end of the spectrum has the overlap of the two contributions and with the underlying 
background. Furthermore, since the splitting between the masses of the /if , fif and that 
of 71 is small, the lower end of the momentum distribution can be as low as 0(1 GeV) 
where the lepton identification efficiency is cut-off by the solenoidal field bending the lepton 
before it reaches the electro-magnetic or the hadron calorimeter [46]. Nevertheless, there 
is sufficient information in this distribution to extract the mass of the 71 particle, using 
the prior information on the /uf and fif masses, obtained by the threshold scan. 

In Fig. |b we show the muon energy distribution after detector simulation. A one 
parameter fit gives an uncertainty on the 71 mass of ±0.19 (stat.) ±0.21 (syst) GeV, where 
the statistical uncertainty is given for 1 ab _1 of data and the systematics reflects the 
effect of the uncertainty on the fi\ masses. The beamstrahlung introduces an additional 
systematics, which depends on the control of the details of the luminosity spectrum. 

4.5 Photon energy spectrum and radiative return to the Z2 

With the e + e~ colliding at a fixed centre-of-mass energy above the pair production thresh- 
old a significant fraction of the KK muon production will proceed through radiative return. 
Since this is mediated by s-channel narrow resonances, a sharp peak in the photon energy 
spectrum appears whenever one of the mediating s-channel particles is on-shell. In case of 
supersymmetry, only Z and 7 particles can mediate smuon pair production and neither of 
them can be close to being on-shell. On the contrary, an interesting feature of the UED 
scenario is that [i\ production can be mediated by Z n and 7 n KK excitations (for n even) as 
shown in Fig. |l]b. Among these additional contributions, the Z2 and 72 exchange diagrams 
are the most important. Since the decay Z2 — * fiifJ-i is allowed by phase space, there will 
be a sharp peak in the photon spectrum, due to a radiative return to the Z2. The photon 
peak is at 

1 f Ml \ 

E, = -E CM y--^)- ( 4 - 5 ) 

On the other hand, M 72 < 2M„ 15 so that the decay 72 — * (J,ifAi is closed, and therefore 
there is no radiative return to 72. Notice that the level 2 Weinberg angle is very small [12] 
and therefore Z2 is mostly VF^-like and couples predominantly to //f and not /if. 

The photon energy spectrum in e + e~ — > fj^fx^^y for R~ l = 1350 GeV, AR = 20 and 
Ecm = 3 TeV is shown in Fig. ||. On the left we show the ISR-corrected theoretical 
prediction from CompHEP while the result on the right in addition includes detector and 
beam effects. It is clear that the peak cannot be missed. 

5. Prospects for discovery and discrimination in other final states 

Previously in section |3] we considered the fi + /j>~]fi final state resulting from the pair pro- 
duction of level 1 KK muons. However, this is not the only signal which could be expected 
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Figure 6: Photon energy spectrum in e + e — > fxf fx x 7 for R 1 = 1350 GeV, AR = 20 and 
Eqm = 3 TeV before (left) and after (right) detector simulation. The acceptance cuts are E 1 > 10 
GeV and 1 < 7 < 179°. The mass of the Z 2 resonance is 2825 GeV. 



in the case of UED. Due to the relative degeneracy of the KK particles at each level, the 
remaining n = 1 KK modes will be produced as well, and will yield observable signatures. 
In those cases, the discrimination techniques which we discussed earlier can still be applied, 
providing further evidence in favor of one model over the other. In this section we compute 
the cross-sections for some of the other main processes of interest, and discuss how they 
could be analysed. 

5.1 Kaluza-Klein leptons 

We first turn to the discussion of the other KK lepton flavors. The KK r- leptons, 7^ , 
are also produced in s-channel diagrams only, as in Fig. |l|, hence the t^t^ production 
cross-sections are very similar to the LiffJ-i case. The final state will be T + T~]fi, and it 
can be observed in several modes, corresponding to the different options for the r decays. 
However, due to the lower statistics and the inferior jet energy resolution, none of the 
resulting channels can compete with the discriminating power of the [i^fi^tf] final state 
discussed in the previous section. 

The case of KK electrons is more interesting, as it contains a new twist. The production 
of KK electrons can also proceed through the t-channel diagram shown in Fig. [7]c. As a 
result, the production cross-sections for KK electrons can be much higher than for KK 
muons. We illustrate this in Fig. ||, where we show separately the cross-sections for SU(2)\y 
doublets (solid lines) and SU(2)w singlets (dotted lines), as a function of R . (For the 
numerical results throughout section ||, we always fix AR = 20.) At low masses (i.e. low 
R^ 1 ) the efe~[ cross-sections can be up to two orders of magnitude larger, compared to the 
case of fiffi^. Another interesting feature is the resonant enhancement of the cross-section 
for R^ 1 ~ 1450 GeV, which is present in either case (e or fx) for the SU(2)w doublets 
(solid lines), but not the SU(2)w singlets (dotted lines). The feature is due to the on-shell 
production of the level 2 Z2 KK gauge boson, which can then decay into a pair of level 
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Figure 7: The same as Fig. H, but for KK electron production 
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Figure 8: ISR-corrected production cross-sections of level 1 KK leptons (e\ in red, fii in blue) at 
ClAC, as a function of Rr 1 . Solid (dotted) lines correspond to SU(2)w doublets (singlets). 



1 KK leptons (see diagram (b) in Figs. |l| and 0). Since the Weinberg angle at the higher 
(n > 0) KK levels is tiny [12], is predominantly an SU(2)w gauge boson and hence 
does not couple to the SU(2)w singlet fermions, which explains the absence of a similar 
peak in the ef and /if cross-sections 4 . 

Because of the higher production rates, the e + e~lfi event sample will be much larger 
and have better statistics than /i + fi~ Ifi. The e + e~ Ifi final state has been recently advertised 
as a discriminator between UED and supersymmetry in [47]. However, the additional t- 
channel diagram (Fig. ^c) has the effect of not only enhancing the overall cross-section, 



but also distorting the differential angular distributions discussed previously in Section 4.1 
and creating a forward peak, which causes the cases of UED and supersymmetry to look 
very much alike. We show the resulting angular distributions of the final state electrons in 
Fig. ^. For proper comparison, we follow the same procedure as before: we choose the UED 
spectrum for R~ l = 500 GeV, which yields KK electron masses as in Table @. We then 



4 One might have expected a second peak closeby due to 72 resonant production, but in the minimal 
UED model the spectrum is such that the decays of 72 to level 1 fermions are all closed. 
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choose a super symmetric spectrum with selectron mass parameters as in Table g. This 
guarantees matching mass spectra in the two cases (UED and supersymmetry) so that any 
differences in the angular distributions should be attributed to the different spins. 

Unlike Fig. [3|, where the underlying shapes of the angular distributions were very 
distinctive (see eqs. Q4.2| ) and ( |4.3| )), the main effect in Fig. ^is the uniform enhancement of 
the forward scattering cross-section, which tends to wash out the spin correlations exhibited 
in Fig. | 

5.2 Kaluza-Klein quarks 

Level 1 KK quarks will be produced in s-channel via diagrams similar to those exhibited in 
Fig. p]. The corresponding production cross-sections are shown in Fig. 1C, as a function of 
R~ l . We show separately the cases of the SU(2)w doublets uf and (if and the SU(2)\y 
singlets uf and df. In the minimal UED model, the KK fermion doublets are somewhat 
heavier than the KK fermion singlets [12], so naturally, the production cross-sections for 
uf and df cut off at a smaller value of R~ x . Since singlet production is only mediated 
by U(l) hypercharge interactions, the singlet production cross-sections tend to be smaller. 
We notice that uf uf is larger by a factor of 2 2 compared to dfdf, in accordance with the 
usual quark hypercharge assignments. 

The observable signals will be different in the case of SU(2)w doublets and SU(2)yy 
singlets. The singlets, uf and df, decay directly to the LKP 71, and the corresponding 
signature will be 2 jets and missing energy. The jet angular distribution will again be 
indicative of the KK quark spin, and can be used to discriminate against (right-handed) 
squark production in supersymmetry, following the procedure outlined in section |4.1| . The 
jet energy distribution will again exhibit endpoints, which will in principle allow for the 



mass measurements discussed in section 4.4. A threshold scan of the cross-section will pro- 
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Figure 10: ISR-corrected production cross-sections of level 1 KK quarks at CLIC, as a function 
ofR- 1 . 



vide further evidence of the particle spins (see section [4^ ). The only major difference with 
respect to the [i^ ~$ final state discussed in section ||, is the absence of the monochro- 



matic photon signal from section [4.5| , since Zi is too light to decay to KK quarks. In 
spite of the many similarities to the dimuon final state considered in section ||, notice that 
jet angular and energy measurements are not as clean and therefore the lepton (muon or 
electron) final states would still provide the most convincing evidence for discrimination. 

The signatures of the SU(2)\y doublet quarks are richer - both uf and df predomi- 
nantly decay to Z\ and which in turn decay to leptons and the LKP [4] . The analogous 
process in supersymmetry would be left-handed squark production with subsequent decays 
to X2 or X\i w hich in turn decay to II and Xv I n principle, the spin information will still 
be encoded in the angular distributions of the final state particles. However, the analysis 
is much more involved, due to the complexity of the signature, and possibly the additional 
missing energy from any neutrinos. 

5.3 Kaluza-Klein gauge bosons 

The ISR-corrected production cross-sections for level 1 electroweak 5 KK gauge bosons 
(Wy , Z\ and 71) at a 3 TeV e + e~ collider are shown in Fig. |ll|, as a function of R" 1 . 
The three relevant processes are W^W^, Z\Z\ and Z\^\ (7171 is unobservable) . In each 
case, the production can be mediated by a t-channel exchange of a level 1 KK lepton, 
while for W+I^ - there are additional s-channel diagrams with 7, Z, 72 and Z2. Z\ and 
W-y are almost degenerate [12], thus their cross-sections cut off at around the same point. 
The analogous processes in supersymmetry would be the pair production of gaugino-like 
charginos and neutralinos. The final states will always involve leptons and missing energy, 
since W x and Z\ do not decay to KK quarks. 

5 The level 1 KK gluon, of course, has no tree-level couplings to e + e~ . 
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In conclusion of this section, for completeness we also discuss the possibility of ob- 
serving the higher level KK particles and in particular those at level 2. For small enough 
R , level 2 KK modes are kinematically accessible at CLIC. Once produced, they will in 
general decay to level 1 particles and thus contribute to the inclusive production of level 1 
KK modes. Uncovering the presence of the level 2 signal in that case seems challenging, 
but not impossible. 

We choose to concentrate on the case of the level 2 KK gauge bosons (V2), which 
are somewhat special in the sense that they can decay directly to SM fermions through 
KK number violating interactions. Thus they can be easily observed as dijet or dilepton 
resonances. In principle, there are two types of production mechanisms for level 2 gauge 
bosons. The first is single production e + e~ — ► V2, which can only proceed through KK 
number violating (loop suppressed) couplings. The second mechanism is e + e~ — > V2V2 pair 
production which is predominantly due to KK number conserving (tree-level) couplings. 



In Fig. |12| we show the corresponding cross-sections for the case of the neutral level 2 gauge 
bosons, as a function of -R -1 , For low values of R" 1 , pair production dominates, but as 
the level 2 gauge boson masses increase and approach Ecm, single production becomes 
resonantly enhanced. Thus the first indication of the presence of the level 2 particles may 
come from pair production events, but once the mass of the dijet or dilepton resonance 
is known, the collider energy can be tuned to enhance the cross-section and study the V2 
resonance properties in great detail. 

6. Conclusions 

Supersymmetry and Universal Extra Dimensions are two appealing examples of new physics 
at the TeV scale, as they address some of the theoretical puzzles of the SM. They also 
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Figure 12: ISR-corrected production cross-sections of level 2 KK gauge bosons at CLIC, as a 
function of 



provide a dark matter candidate which, for properly chosen theory parameters, is consistent 
with present cosmology data. Both theories predict a host of new particles, partners of 
the known SM particles. If either one is realised in nature, the LHC is expected to observe 
signals of these new particles. However, in order to clearly identify the nature of the new 
physics, one may need to contrast the UED and supersymmetric hypotheses at a multi-TeV 
e + e" linear collider such as CLIC 6 . In this paper we studied in detail the process of pair 
production of muon partners in the two theories, KK-muons and smuons respectively. We 
used the polar production angle to distinguish the nature of the particle partners, based 
on their spin. The same analysis could be applied for the case of other KK fermions, as 
discussed in section ||. 

We have also studied the accuracy of CLIC in determining the masses of the new 
particles involved both through the study of the energy distribution of final state muons and 
threshold scans. An accuracy of better than 0.1% can be obtained with 1 ab" 1 of integrated 
luminosity. Once the masses of the partners are known, the measurement of the total cross 
section serves as an additional cross-check on the hypothesized spin and couplings of the 
new particles. A peculiar feature of UED, which is not present in supersymmetry, is the 
sharp peak in the ISR photon energy spectrum due to a radiative return to the KK partner 
of the Z. 

The clean final states and the control over the centre-of-mass energy at the CLIC 
multi-TeV collider allows one to unambiguously identify the nature of the new physics 
signals which might be emerging at the LHC already by the end of this decade. 

6 Similar studies can also be done at the ILC provided the level 1 KK particles are within its kinematic 
reach. Since precision data tends to indicate the bound > 250 GeV for the case of 1 extra dimension, 
one would need an ILC center-of-mass energy above 500 GeV in order to pair-produce the lowest lying KK 
states of the minimal UED model. 
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